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The coarsening kinetics for the beading instability for liquid contained in a groove with convex 
curved sides (for example between a pair of parallel touching cylinders) is considered as an open 
channel flow problem. In contrast to a V-shaped wedge or U-shaped microchannel, it is argued 
that droplet coarsening takes place by viscous hydrodynamic transport through a stable column of 
liquid that coexists with the droplets in the groove at a slightly positive Laplace pressure. With 
some simplifying assumptions, this leads to a tV' growth law for the characteristic droplet size as 
a function of time, and a t~ i ^ law for the decrease in the droplet line density. Some remarks are 
also made on the spreading kinetics of an isolated drop deposited in such a groove. 

PACS numbers: 47.55.nb, 47.20.Dr 


I. INTRODUCTION 

Open channel flow problems have attracted much in¬ 
terest not only because of possible applications in mi¬ 
crofluidics m , but also because of their relevance to 
liquids spreading on topographically patterned surfaces 
such as human skin [TUUT2] . The paradigmatic case of 
spreading in a V-shaped wedge has been analysed both 
when liquid is supplied by a reservoir uni na, and in 
the starved (no reservoir) situation [14] . Various aspects 
of these predictions have been confirmed experimentally 
EDUE]. Flows in U-shaped channels ( i. e. with concave 
sides) have also been considered since such channels are 
easily micro-machined and are thus relevant for microflu¬ 
idics applications mini- In the present study, I revisit 
the problem, in the context of a groove with convex sides, 
such as that formed between a pair of parallel touching 
cylinders (Fig. [IJ. In addition to possible microfluidics 
applications, interest in this problem is motivated by con¬ 
sideration of oily soil spreading along the fibres in tex¬ 
tile yarns in woven fabrics m ■ It also has other poten¬ 
tial technological relevance, for example to molten solder 
wicking in stranded copper wires and braids. 

This problem throws up some interesting aspects not 
found in the previous cases. A groove with convex curved 
sides supports a stable uniform liquid column at low 
loads, but at higher loading a uniform liquid column may 
display a beading instability in which it breaks up into a 
string of droplets, similar to Rayleigh’s observation of the 
breakup of a thread of treacle on a paper surface fIB] . 
Although the statics are by now quite well understood 
[5] ;21 O [T5] lTj [TSj , I shall argue here that the loading 
duality uniquely differentiates the present case from V- 
shaped wedges and U-shaped microchannels, since the 
droplets (beads) can coarsen by mass transport through 
connecting liquid columns. 

Additionally, the loading duality implies an isolated 
droplet deposited into such a groove may show several 


kinetic spreading stages as it empties into the unfilled 
regions. This aspect will be discussed at the end. 


II. WICKING EQUATION 


To start with, consider the general problem of open 
channel flow in a channel of arbitrary cross section, and 
let A{x , t) be the cross section occupied by liquid. A local 
mass conservation law holds mniiig, 


dA d{Av) 
dt dx 


(1) 


Herein, the mean flow rate v satisfies a Hagen-Poiseuille 
(HP) law, v = — {k/rj) dp/dx, in which k is the per¬ 
meability (a quantity with units of length squared, cf. 
Darcy’s law), 77 is viscosity, and p = p{A) is the (loading- 
dependent) Laplace pressure. Combining the HP law 
with Eq. ([!]) gives what can perhaps be called the wicking 
equation, 


dA d / Ak dp <9A\ 
dt dx V rj dA dx ) ' 

This is the basis for much of the subsequent development, 
and also codifies the statics via stability analysis. It gen¬ 
erally has the character of a non-linear diffusion equa¬ 
tion. I have assumed that the occupied cross section A is 
weakly varying with x, so the contribution to the Laplace 
pressure from the interface curvature in the longitudinal 



FIG. 1. Liquid in a groove between parallel touching cylin- 
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direction can be neglected. Note that the permeability 
also depends on the loading, so that k = k(A). 

Insight can be gained by linearising about the uni¬ 
formly loaded static solution, viz. A = A 0 and v = 0. Let 
us write A/A 0 = 1 + e(x,t)- Then de/dt = D e g d 2 e/dx 2 
where D e g = (Ak/rj) dp/dA is an effective diffusion co¬ 
efficient, evaluated at A = A 0 . If this is positive, then 
perturbations will die away and the liquid will be self¬ 
levelling. If it is negative, then perturbations will grow 
indicative of the aforementioned beading instability. It 
is clear that the behavior depends, not on the sign of 
Laplace pressure p , but rather on the sign of dp/dA. 
Physically, if dp/dA > 0, an overfilled region will have a 
higher Laplace pressure than an underfilled region, and 
the liquid will flow to even things out. On the other hand, 
if dp/dA < 0, liquid will flow from underfilled regions into 
overfilled regions, magnifying the initial imbalance. The 
interesting and unusual property of a groove with convex 
curved sides is that both situations occur, depending on 
the loading. 

III. STATICS 

A specific example of the loading duality is provided by 
the groove between a pair of parallel touching cylinders, 
shown in Figs. |T] and [ 2 j Since I neglect the interface cur¬ 
vature in the longitudinal direction, the transverse profile 
of the free surface is characterised by an arc of a circle 
with radius R. Taking R > 0 to indicate the surface is 
convex outwards, one has p = 7 /R where 7 is surface ten¬ 
sion. In this problem it is convenient [5] to parametrise 
the loading by the wrapping angle a, shown in Fig. |Tj 
Note that the area is a monotonically increasing function 
of a , so dA/da > 0. 

Elementary trigonometric arguments, first presented 
to my knowledge by Princen he show that 

a, = cos( 6 > + a) 

R 1 — cos a 

where 6 is the liquid-solid contact angle. This result holds 
for both convex and concave interfaces, and the sign has 
been inserted in accord with the above convention that 
R is positive if the interface is convex. 

From this we find (after a little rearrangement) 

dR _ R 2 cos(^a + 9) 
da 2 a sin 3 (I a) 

Since p oc a/R, Eqs. ([3]) and Q show that p increases 
through zero at a + 6 = | 7 r to reach a weak maximum 
at ^a + 9 = I- 7 T (i. e. a = ir — 29). A specific example 
is shown in Fig. Sb for a contact angle 9 = 30°. The 
zero crossing is at a = 60° and the Laplace pressure 
maximum is at a = 120°. Fig. I3> shows a selection of the 
corresponding filling states. The stability requirement 
that dp/dA oc dp/da > 0 indicates that filling states 
with a > 120° are unstable with respect to the above- 
mentioned beading instability [ 20 ;. 




FIG. 2. (a) The dimensionless Laplace pressure as a function 
of the filling angle, with liquid-solid contact angle 6 = 30°. (b) 
Filling stages in this problem: the thick solid line is the max¬ 
imum Laplace pressure at a = 120°; uniform loading above 
this (dashed lines) is predicted to be unstable. 


IV. KINETICS 
A. Beading instability 

As we have seen for the case of parallel touching cylin¬ 
ders, the Laplace pressure in a groove with convex sides 
may show a maximum as a function of the loading. Be¬ 
yond the maximum, dp/dA < 0, and therefore a beading 
instability arises in which the uniformly loaded state is 
unstable towards the growth of perturbations. But first, 
what could be the final state in such a situation? In prin¬ 
ciple one can have two different states of loading at same 
Laplace pressure, which can therefore be in coexistence, 
for example points P and Q in in Fig. However the 
higher loaded state is always in the unstable region. The 
logical conclusion is that the excess liquid is expelled into 
a large droplet that sits somewhere on the two cylinders, 
coexisting with a stable column of liquid in the groove. 
In the case of parallel touching cylinders, the wrapping 
angle for this stable liquid column would have to satisfy 
— 9 < a < tt — 29. The first inequality arises because 
the liquid column coexists with a large drop at a (weak) 
positive Laplace pressure [2Tj. The second inequality is 
the column stability condition. 

I turn now to the coarsening kinetics. A uniformly 
overloaded state will certainly break up into a string of 
droplets, but in a groove with convex curved sides these 
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FIG. 3. Mean droplet volume as a function of time. Droplets 
coarsen by viscous hydrodynamic transport along the con¬ 
necting liquid columns (diagrammatic inset), according to the 
rules prescribed in the Appendix. Results from 10 indepen¬ 
dent simulation runs are shown. Each simulation was ini¬ 
tialised with 5 x 10 4 equispaced droplets, with random sizes 
taken from a uniform distribution Vo < V < 2Vo- Droplets 
are removed when V < 0.9 Vo, and their liberated contents 
are added to the liquid column (/3 = 1 in the model). The 
time step parameter was t = 0.02. 


droplets are always connected by liquid columns with 
a finite filling depth according to the above argument. 
Therefore the larger droplets can eat the smaller ones, by 
transporting liquid along the connecting liquid columns. 
This stands in contrast to the V-shaped wedge for exam¬ 
ple, where spatially separated droplets are disconnected 
m HH123] and the droplet population has to coarsen by 
some other mechanism, such as a prewetting film EU- 
[26] , or transport through the vapor phase in the case of 
a volatile liquid. 

With a simplifying assumption about how the Laplace 
pressure depends on the droplet size, aspects of the 
Lifshitz-Slyozov-Wagner (LSW) theory of droplet coars¬ 
ening m can be adapted to the present case. The sim¬ 
plifying assumption is that the Laplace pressure scales in¬ 
versely with droplet size as p oc V -1//3 , where V = R 3 is 
the droplet volume. This is likely to be true only asymp¬ 
totically [251 125]. Nonetheless, let us suppose that this 
is true, and the situation has evolved so that there is 
a string of droplets sitting in the groove. The follow¬ 
ing mean-field scaling ansatz predicts how the character¬ 
istic droplet size R, and characteristic spacing between 
droplets L, evolve with time. 

First, the mass flux between adjacent droplets will be 


J 


a 4 A P 
V L 


(5) 


This just expresses the HP law in scaling form. Shown 
here is the prefactor for the case of parallel touching 
cylinders, where the fourth power of the cylinder radius a 
arises from the product of the Darcy permeability k ~ a 2 
and the cross sectional area of the connecting liquid col¬ 
umn A ~ a 2 . The analysis does not depend on this 


FIG. 4. Scale invariance of droplet size distribution. His¬ 
tograms are computed when there are 10 000, 5000, 1000, 
500, 300 and 100 droplets remaining out of an initial 5 x 10 4 , 
combining data from 100 independent simulation runs. Sim¬ 
ulation parameters as for Fig. [3] 


particular geometry though. 

The mass flux is incorporated into a local mass conser¬ 
vation law, cf. Eq. (jTJ) , which governs the growth of the 
mean droplet volume, 


dV 

dt 


~ J. 


( 6 ) 


At the same time there is a global mass conservation law 
which relates the mean droplet size to the mean spacing, 


V 

— ~ oj ~ const, 


(7) 


where w is the mass per unit length [3D]. 

Combining Eqs. § 0' together with V = R 3 and 
A p ~ j/R, shows that 


R 2 


dR 

dt 


7 a 4 uj 
r]R 4 


This integrates to 


^7 a 3 wt^ V 7 


( 8 ) 

(9) 


Thus the prediction is that the mean droplet size should 
grow as f 1 / 7 and the droplet line density (i. e. 1/L ~ R~ 3 ) 
should diminish as f -3 / 7 . 

The mean-field assumption is questionable given the 
one-dimensional nature of the coarsening problem. To 
investigate specifically just this aspect, I undertook nu¬ 
merical simulations using the model described in the Ap¬ 
pendix. This confirms that the mean-field scaling ansatz 
does indeed predict how the mean droplet size grows with 
time (Fig. [3]), and also demonstrates LSW-like scale in¬ 
variance for the droplet size distribution (Fig. |4|) . 

The origin of the mean-field behaviour is apparent if 
one examines the equal-time droplet size correlation func¬ 
tion, C{\i — j |), shown in Fig. [ 5 ] As can be seen the 
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FIG. 5. Equal-time size correlation function, C(\i — j\), com¬ 
puted for the same set of simulations used for Fig. [4] The 
inset shows the time dependence of the depth of the nearest- 
neighbour minimum. Error bars are from block averaging (10 
blocks x 10 runs). 

only significant correlation appears to be between near¬ 
est neighbours, at \i—j\ = 1. This indicates that droplets 
which are larger than average tend to be adjacent to 
droplets which are smaller than average, but apart from 
this no significant long-range correlations develop. 

To see the origin of this nearest-neighbour correla¬ 
tion hole, consider an artificial situation in which one 
large droplet sits in a uniform string of equisized smaller 
droplets. Away from the large droplet, the ambient 
Laplace pressure is uniform and no coarsening takes 
place. However the large droplet has a sub-ambient 
Laplace pressure and so starts to draw in material from 
its immediate neighbours. This causes the immediate 
neighbours to shrink, and increases their Laplace pres¬ 
sure relative to the ambient background. The next- 
nearest neighbours then see that these droplets have 
started to shrink, and so they in turn start to grow. The 
process continues, and it is easy to see that it will gener¬ 
ate a staggered array of droplet sizes, with a concomitant 
nearest-neighbour negative size correlation. 

The nearest-neighbour correlation hole diminishes 
somewhat as time progresses, but eventually appears to 
settle down to a value C( 1) ~ —0.15 (Fig. [5] inset). One 
should point out that the coarsening dynamics in the 
model is quite subtle, for instance, the largest droplet at 
time t may not necessarily be the largest droplet at some 
later time t' > t. This is because the growth rate of a 
droplet depends not only on its size but also on the sizes 
of its neighbours, and how far away they are. 


B. Droplet spreading 

The remaining point of kinetic interest concerns the 
fate of a droplet of liquid deposited onto an initially 
empty groove. The droplet will start to empty into 
the groove, presumably driving a Bell-Cameron-Lucas- 


Washburn (BLCS) type flow from what is in effect a 
shrinking droplet reservoir [Iff, :T5j [19]. This should per¬ 
sist all the way until the loading falls below the Laplace 
pressure maximum. Past this point, by analogy to the 
V-shaped wedge na, one expects the spreading rate to 
slow down since the reservoir has been exhausted. 

For the case of parallel touching cylinders, a second 
power law appears at a very late stage where every¬ 
where a <C 1. In this limit, the wetted portion of the 
groove has shrunk to a narrow fissure with a width of 
the order R ~ a 2 a (see Eq. ([ 3 ]) in the limit a —> 0) 
and a depth of the order aa. One therefore expects 
A ~ a 3 a 2 , and presumably k ~ a 4 a 2 since the permeabil¬ 
ity should largely be determined by the width. Similar 
to the coarsening kinetics problem, a scaling analogue 
of the HP law can be introduced. In the present case 
this is dL/dt ~ (k/rj) x Ap/L where L is the length 
of the wetted portion of the groove and A p ~ 7 /R is 
the Laplace pressure. Substituting the above scaling ex¬ 
pressions gives dL/dt ~ 'ya 2 a/( , qL). An additional con¬ 
straint comes from the analogue of the global mass con¬ 
servation law in Eq. 0, namely that the total droplet 
volume fl ~ a 3 a 2 L should be conserved. Eliminating a 
between this volume constraint and the HP scaling law 
yields dL/dt ~ 7 H 2 / 3 /(p^^L 5 / 3 ). This integrates to 
the final rather esoteric result 

L ~ (7t/r?) 3/8 H 1/4 a- 1/8 . (10) 

In other words, the initial L ~ t 1 / 2 spreading law (BLCS) 
should weaken when the droplet reservoir vanishes, and 
eventually enter an L ~ f 3 / 8 power law in the final 
starved state. 


V. DISCUSSION 

I have argued that the kinetic aspects of spreading and 
(de)wetting for a liquid contained in a groove with con¬ 
vex curved sides presents some unique aspects when com¬ 
pared, for example, to a V-shaped wedge or a U-shaped 
microchannel. The novel aspects arise from an underly¬ 
ing loading duality, wherein a liquid column is stable at 
low loading, but becomes unstable at higher loading. 

For the case of uniform loading above the critical load¬ 
ing threshold, a beading instability should be observed 
in which the liquid column breaks up into a string of 
droplets, which subsequently coarsen by mass transport 
along connecting liquid columns. With a simplifying as¬ 
sumption about the dependence of the Laplace pressure 
on the droplet volume, a mean-field scaling ansatz in¬ 
dicates that the droplet size and line density scale with 
non-trivial power laws in time, and simulations show that 
this is not destroyed by the one-dimensional nature of the 
problem. 

Of course, the simplifying assumption about the 
Laplace pressure scaling does not hold in reality, and in 
general I would expect the clean power law behaviour to 
be modified by finite-droplet-size effects, which may be 
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quite persistent. Some general predictions should be ro¬ 
bust however, such as the relatively slow droplet growth 
via transport along connecting liquid columns, and the 
negative correlation between nearest neighbour droplet 
sizes shown in Fig. [5j These could perhaps be tested in 
an electrowetting experiment min- I should caution that 
the predicted slow coarsening kinetics may be overtaken 
by other, ultimately faster, mechanisms. For example 
transport through the vapour phase for a volatile liquid 
may ultimately lead to a i? ~ t 1 / 3 coarsening law, as in 
LSW theory [127.. 

Another prediction arising from the loading duality is 
that a droplet deposited in the groove should show a 
staged spreading kinetics, starting with the classic Bell- 
Cameron-Lucas-Washburn law as the droplet initially 
acts like a reservoir, slowing when the loading falls every¬ 
where below the critical value, and possibly ending with 
a new power law in the final starved state. Again, this 
may perhaps be probed experimentally. 


Appendix A: Simulation of coarsening kinetics 


I introduce a simplified model of the coarsening kinet¬ 
ics to test specifically the mean-field ansatz presented 
in the main text. In the model, I consider a one- 
dimensional string of i = 1... IV droplets, connected by 
liquid columns as indicated in the lower inset in Fig. [3] 
Big droplets grow, and small droplets shrink, under the 
influence of viscous hydrodynamic transport through the 
liquid columns. This is driven by differences in the 
Laplace pressure between neighbouring droplets. To es¬ 
tablish a system of kinetic equations for the droplet sizes, 
I suppose that the Laplace pressure in the i-th droplet 
is proportional to Vi , where V) is the droplet volume 
(this is the simplifying assumption mentioned in the main 
text). According to the HP law, the Laplace pressure dif¬ 
ference drives a mass flux through the connecting liquid 
column as, cf. Eq. |5]), 


Ji = 


y- 1,3 -v- + T 

Li 


(Al) 


In this Li is the distance between the i-th and (i + l)-th 
droplets, and all other material properties in the problem 
have been subsumed into the definitions of length, volume 
and time. Note that V)+i > V) implies Ji > 0, so that liq¬ 
uid flows from smaller droplets to larger droplets. Given 
the fluxes, mass conservation dictates that, cf. Eq. ([6]), 

^ = Ji~\ - J t . (A2) 


Eqs. (Al) and (A2) are the required set of non¬ 
linear kinetic equations. Since they predict that droplets 
shrink, as well as grow, we need a rule which governs how 
shrinking droplets can disappear. At this point it is con¬ 
venient to introduce a fiducial volume Vo ~ Rq, where 
the fiducial length R 0 is set by the height of the con¬ 
necting liquid column (in the simulations, Vo = i? 0 = 1). 
A simple rule for shrinking droplets is that they vanish 
when V. < aVg. If this happens, the droplet is removed 
and the distance between the remaining droplets is set 
equal to + Li + f3V i ' , where the third term is an 
ad hoc correction for the length contributed by the van¬ 
ished droplet (taking /3 as a free parameter). 


As an initial condition I set Vi = ?yVo where r t is a 
random number chosen from a uniform distribution, 1 < 
r,; < r m . The droplets are initially equispaced, with Li = 
Rq. Periodic boundary conditions are imposed. 


The droplet vo lumes are evolved according to 
Eqs. ( |A1[ ) and (A2), alongside the above rule for remov¬ 
ing droplets which become too small. Eqs. (Al) and (A2) 
are integrated using a simple, adaptive, Euler-type for¬ 
ward finite difference scheme, with a time step At chosen 
such dVi/dt x At/Vi < e, in other words so that the frac¬ 
tional change in any droplet size does not exceed e in any 
time step. For the reported simulations I used a = 0.9 
and /3 = 1 for the vanishing rule, r m = 2 for the max¬ 
imum initial drop size relative to Vo, and e = 0.02 for 
the choice of time step. I have checked the results are 
insensitive to these choices. 


As time evolves, the larger droplets grow at the ex¬ 
pense of the smaller droplets, and the smallest droplets 
shrink and vanish. Eventually the simulation stops when 
there is one large final droplet (N = 1). I monitor the 
mean droplet volume (V) = (l/IV) YliLi Vi as a func¬ 
tion of time (note that N changes as droplets disappear), 
and at selected time points record the drop size distri¬ 
bution. Typical results, aggregated from multiple sim¬ 
ulation runs from independent starting points, are sum¬ 
marised in Figs. [3] and [4] 

I also calculate periodically the equal-time correlation 
function 


|) 


(AVi AVj) 

(AH 2 ) 


(A3) 


where AV) = V) — (V) is the deviation from the mean 
droplet volume, at time t. This function is shown in 
Fig.m evaluated at various points in the simulation. 
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